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Recent progress in computing four-loop massive correlators * 
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We report about recent progress in computing four- loop massive correlators. The expansion of these correlators 
in the external momentum leads to vacuum integrals. The calculation of these vacuum integrals can be used to 
determine Taylor expansion coefficients of the vacuum polarization function and decoupling functions in pertur- 
bative Quantum chromodynamics. New results at four-loop order for the lowest Taylor expansion coefficient of 
the vacuum polarization function and for the decoupling relation are presented. 



1. Introduction 

Two-point correlators have been studied in 
great detail in the framework of perturbative 
quantum field theory. Due to simple kinematics 
(only one external momentum) even multi-loop 
calculations can be performed. The results for all 
physically interesting diagonal and non-diagonal 
correlators and including full quark mass depen- 
dence are available up to O (a%) [1-3]. 
At four-loop order the two-point correlators can 
be considered in two limits. In the high en- 
ergy limit massless propagators need to be cal- 
culated and in the low energy limit vacuum dia- 
grams (tadpole integrals without dependence on 
the external momentum) arise. The evaluation of 
these massive tadpoles in three-loop approxima- 
tion has been pioneered in ref. [4] and automated 
in ref. [5] . 

Similar to the three-loop case, the analytical eval- 
uation of four-loop tadpole integrals is based 
on the traditional Integration-By-Parts (IBP) 
method. In contrast to the three-loop case the 
manual construction of algorithms to reduce ar- 
bitrary diagrams to a small set of master integrals 
is replaced by Laporta's algorithm [6,7]. In this 
context the IBP identities are generated with nu- 
merical values for the powers of the propagators 
and the irreducible scalar products. In the next 
step, the resulting system of linear equations is 



then solved in the next step by expressing system- 
atically complicated integrals in terms of simpler 
ones. The resulting solutions are then substituted 
into all the other equations. 

This reduction has been implemented in an au- 
tomated F0RM3 [8, 9] based program in which 
partially ideas described in ref. [7, 10, 11] have 
been implemented. The rational functions in the 
space-time dimension d, which arise in this proce- 
dure, are simplified with the program FERMAT [12]. 
The automated exploitation of all symmetries of 
the diagrams by reshuffling the powers of the 
propagators of a given topology in a unique way 
strongly reduces the number of equations which 
need to be solved. 

In general, the tadpole diagrams encountered 
during our calculation contain both massive and 
massless lines. In contrast, the computation 
of the four-loop /3-functions can be reduced to 
the evaluation of four-loop tadpoles composed of 
completely massive propagators. These special 
cases have been considered in [11, 13, 14]. 

The outline of this paper is as follows. In sec- 
tion |21 we discuss the calculation of the lowest 
expansion coefficient of the vacuum polarization 
function and present the results at four-loop or- 
der using methods as described above. In section 
01 we present new results for the decoupling rela- 
tion at four- loop order in perturbative QCD. Our 
conclusions are presented in section 01 



"Talk given at the 12th International Conference on Quan- 
tum Chromodynamics, Montpellier, 4-8th July 2005 
tOn leave from Institute for Nuclear Research of the Rus- 
sian Academy of Sciences, Moscow, 117312, Russia. 



2. Vacuum polarization function 

The vacuum polarization tensor Il tll/ (q) is de- 
fined as 



IL^(q)=i dxe i<lx (0\Tf(x)f(0)\0), 



(1) 



where g M is the external momentum and j M is 
the electromagnetic current of a heavy quark with 
mass to /j. The tensor Ii tiv {q) can be expressed by 
a scalar function, the vacuum polarization func- 
tion H(q 2 ) through 



n**" («) = («"«" 



<?vi 



U{q 2 )+q»q»n L {q 2 ).{2) 



The longitudinal part IIl(<7 2 ) vanishes due to the 
Ward identity. The polarization function II(g 2 ) is 
related to the experimentally measurable _R-ratio 
R(s) through the dispersion relation: 



n(<7 2 ) = n( g 2 = o) 



12 7T 2 



ds- 



R(s) 



s (s — q 2 ) 



(3) 



Performing the n-th derivative of eq. (|3J) with re- 
spect to q 2 at q 2 = one obtains the moments 
M. n , which can be determined experimentally: 



M 



exp 






127T 2 / d . 



The derivatives of the polarization function on 
the rhs are related to the Taylor expansion coef- 
ficients C n : 



n( 9 2 ) 



3Q| 

16 7T 2 



Z^Cn 



(4) 



J!>0 



(z = q 2 /(4rn 2 l )) which can be calculated in per- 
turbative QCD. The first and higher derivatives 
are important for a precise determination of the 
charm- and bottom-quark mass (see e.g. [15]). 
But also the lowest expansion coefficient Co has 
an interesting physical meaning: it relates the 
coupling of electromagnetic interaction in differ- 
ent renormalization schemes. In the case of QED- 
on-shell-renormalization the residue of the pho- 
ton propagator is one and the electrical charge 
e coincides with the classical limit. If one per- 
forms renormalization in the MS-scheme one ob- 
tains a relation between the coupling constant of 
the electromagnetic interaction o cm = e 2 /{An) in 



QED-on-shcll-renormalization and the coupling 
constant a em = e 2 /(47r) in the MS-scheme: 



a em 



l + e 2 II(<7 2 = 0) 



(5) 



For massive quarks, interacting through gluons, 
II(<7 2 = 0) has been evaluated in ref. [1]. At three- 
loop order in perturbative QCD this relation has 
already been determined in ref. [1] . For the QED 
case the corresponding result was calculated in 
ref. [4]. _ 

The first Taylor coefficient Co has been calcu- 
lated using the methods described in section ^ 
All tadpole diagrams were expressed through the 
set of 13 master integrals shown in figure ^ 
These master integrals have been calculated in 




Figure 1. Master integrals. The solid lines de- 
note massive lines, whereas the dashed lines de- 
note massless lines. 



refs. [16-18]. Inserting the master integrals into 
the lowest Taylor coefficient of the polarization 
function and performing the renormalization in 



7 3- 1 



the MS-scheme one obtains the following result 
for the expansion coefficient Co(/i = fnu)' 

Co = (—) 1.4444+ (—} (0.3714 nj + 1.725) 

+ (— ) (0.0257raf- 1.186 nj- 1.955), (6) 

where a s is the strong coupling constant in the 
MS-scheme and the symbol m denotes the num- 
ber of light quarks considered as massless. The 
part of the four-loop contribution proportional to 
the number nl of active quarks has been calcu- 
lated previously up to the first physical moment 
C\ in ref. [16]. The terms proportional to a- 
are even known to all orders j [19]. 
In ref. [20, 21] the terms being proportional to 
nf have been calculated through an independent 
numerical method, which is based on the Pade 
approximation [22-24] . These results are in com- 
plete agreement with those of eq. © . 

3. Decoupling relation 

The masses of the known quark species differ 
vastly in their magnitude. Often the mass of a 
heavy quark h is much larger than the character- 
istic momentum scale y/s of the reaction under 
consideration. In such a two interrelated prob- 
lems appear when using an MS-like renormaliza- 
tion scheme. 

First, two large but quite different mass scales, 
y/s and nih, lead to two different types of poten- 
tially dangerously large logarithms of \[sj [i an d 
mh/n and the standard trick of a proper choice 
of the renormalization scale \i is no longer appli- 
cable; 

Second, according to the Appclquist-Carazzonc 
theorem [25] the effects due to heavy particles 
eventually should in general 'decouple' from low- 
energy physics. However, a peculiarity of mass- 
independent renormalization schemes is that the 
decoupling theorem does not hold in its naive 
form for theories renormalized in such schemes. 
The effective QCD action that appears will not 
be not canonically normalized. Large mass log- 
arithms appear, when one calculates a physical 
observable. 
Fortunately, both problems are controlled once 



the expansion parameters are properly choosen 
and renormalization group improvement is per- 
formed [26-29]. To be specific, consider QCD 
with ni — rif — 1 light quarks and one heavy quark 
h with mass to;,. The effective coupling constant 
a' s is then expressed in terms of the one of the 
full theory through: 

«'»(m) = a s{p) Cfl («»(/■*)) x), x = log(n 2 /m 2 h ) , (7) 

where £ g is the decoupling function and m^di) 
is the MS running mass of the heavy quark. The 
decoupling function is known up to three-loop or- 
der [30,31]. Its calculation can be reduced to 
the solution of vacuum integrals [30]. Thus the 
methods described in section Q can be applied to 
calculate the decoupling function at four-loop or- 
der. The same master integrals as those shown 
in figure ^ appear. For a renormalization scale 
/i = rrih(fJ-) one obtains the following numerical 
result 3 : 



Q = 1 + ( — J 0.1528 

(rv \ 3 
—J (0.9721- 0.0847 m) 



+ 



(8) 



(— ) (5.1703- 1.001 m- 0.0220 nf) 



The knowledge of this decoupling function is of 
phenomenological importance, because it allows 
the determination of a s (Mz) at the Z-Boson 
scale through evolution of the measured value of 
a s (m T ) at the r-lepton scale. A careful analysis 
of the effects of four-loop running and three-loop 
matching on the extraction of a s (Mz) defined 
for 5 active quark flavors from a, (to t ) defined 
for 3 active quark flavors has been recently per- 
formed [34]. We have checked that the inclusion 
of the newly computed four-loop matching condi- 
tion leads to a further reduction of the theoretical 
error from the evolution. 



3 After the completion of our calculation we have been in- 
formed that an independent calculation of the matching 
function £? at the four-loop level has been recently finished 
by Y. Schroder and M. Steinhauser. The same is true for 
Co of eq. JHJ. Their results are in full agreement with 
ours. Both calculations have been made reecntly available 
in [32,33]. 



4. Conclusion 

We have calculated the lowest Taylor expansion 
coefficient of the vacuum polarization function 
and the decoupling relation at four-loop order in 
perturbative QCD. The lowest Taylor expansion 
coefficient relates the electromagnetic coupling 
in the on-shell and in the MS renormalization 
schemes. The decoupling relation is important 
for the determination of a s (Mz) through evo- 
lution from the measured value of a s (m T ). 
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